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ABSTRACT 

When CMB data are used to derive cosmological parameters, their very choice 
does matter: some parameter values can be biased if the parameter space does 
not cover the "true" model. This is a problem, because of the difficulty to 
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parametrize Dark Energy (DE) physics. We test this risk through numerical 
experiments. We create artificial data for dynamical or coupled DE models and 
then use MCMC techniques to recover model parameters, by assuming a constant 
DE state parameter w and no DM-DE coupling. For the DE potential considered, 
no serious bias arises when coupling is absent. On the contrary, uJo^c-, and thence 
Ho and f2o,m, suffer a serious bias when the "true" cosmology includes even just 
a mild DM-DE coupling. Until the dark components keep an unknown nature, 
therefore, it can be important to allow for a degree of freedom accounting for DM- 
DE coupling, even more than increasing the number of parameters accounting 
for the w{a) behavior. 

Subject headings: Dark Matter - Dark Energy - Cosmic Microwave Background 
- Cosmological Parameters 

1. Introduction 

Scarce doubts remain that Dark Energy (DE) exists. Not only SNIa data indicate an 
accelerated cosmic expansion (Perlmutter et al. 1997, 1998, Riess et al. 1998, Foley et al. 
2007); also CMB and deep sample data show a clear discrepancy between the total density 
parameter approaching unity, and the matter density parameter Vto-m ~ 0.25-0.3 (see, 
e.g., Spergel et al. 2007). DE covers this gap; its state parameter w = pde/ Pde must approach 
— 1 today, so apparently excluding that DE is made of free particles (po.de, Po,de'- DE pressure, 
energy density). The true nature of DE is however still elusive; a false vacuum and a self- 
interacting scalar field are among the most popular hypotheses for it (Wetterich 1988, Ratra 
& Peebles 1988). 

In this paper we explore some possible consequences of our poor knowledge of DE nature. 
In particular we test the risk that other cosmological parameter estimates are biased by a 
inadequate parametrization of the DE component. We shall see that this risk is real. 

Theoretical predictions had an astonishing success in fitting CMB data. For instance, 
the SW effect, predicting low-Z Ci data, or primeval compression waves, predicting Ci peaks 
and deeps, were clearly detected. There is little doubt that we are exploring the right range 
of models. 

When we investigate DE nature through CMB data, we must bear in mind that they 
were mostly fixed at a redshift z ~ 1100 , when the very DE density should be neghgible, 
and so affects peak and deep positions indirectly, through the values of uJo,c = ^o,ch^ and 
cUofi = ^o,bh^ (here — 100 /ikm/s/Mpc is the present Hubble parameter; Qq, ^o,c, ^o,b 
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are the present total, CDM, baryon density parameters). Later information on DE state 
equation, conveyed by the ISW effect, is seriously affected by cosmic variance and often relies 
on the assumption that a single opacity parameter r can account for reionization, assumed 
to be (almost) instantaneous. Accordingly, if we assume dynamical DE (DDE), due to a 
scalar field self interacting through a potential V {(/)), CMB data allow to exclude some 
interaction shape, e.g. Ratra-Peebles (1988) potentials with significantly large A energy 
scales, but hardly convey much information on potential parameters. 

In spite of that, when we choose a DE potential or a specific scale dependence of the DE 
state parameter w(a), wc risk to bias the values of other cosmological parameters, sometimes 
leading to premature physical conclusions. An example is the value of the primeval spectral 
index for scalar fluctuation Ug. Using WMAP3 data (Spergel et al. 2007) and assuming 
a ACDM cosmology, the value = 1 is "excluded" at the 2-a confidence level. On the 
contrary, Colombo & Gervasi (2007) showed that this is no longer true in a DDE model 
based on a SUGRA potential (Brax & Martin 1999, 2001; Brax, Martin & Riazuelo 2000), 
whose hkelihood was the same of ACDM. 

The risk that our poor knowledge of DE nature biases parameter determination is even 
more serious if DM-DE couphng is allowed. Coupled DE (CDE) cosmologies were studied 
by various authors (see, e.g., Wetterich 1995, Amendola 2000, Bento, Bertolami & Sen 2002, 
Maccio et al. 2004). 

While DDE was introduced in the attempt to ease DE fine-tuning problems, CDE tries 
to ease the coincidence problem. Let us then parametrise the strength of DM-DE coupling 
through a parameter j3, defined below. When f3 is large enough, DM and DE scale (quasi) 
in parallel since a fairly high rcdshift. In turn this modifies the rate of cosmic expansion 
whenever DM and/or DE contributions to the total energy density are non-negligible, so 
that limits on /? can be set through data. 

The range allowed {f3 < 0.10-0.12; Mainini, Colombo & Bonometto 2005, Mainini 
& Bonometto 2007, see also Majerotto, Sapone & Amendola 2004, Amendola, Campos & 
Rosenfeld, 2006), unfortunately, is so limited that DM and DE are doomed to scale differently, 
but in a short redshift interval. Clearly, this spoils the initial motivation of coupling, but, 
once the genie is outside the lamp, it is hard to put him back inside: even though the coupling 
solves little conceptual problems, we should verify that no bias arises on the other parameters, 
for the neglect of /3's consistent with data. This is far from being just a theoretical loophole, 
the still unknown physics of the dark components could really imply the presence of a mild 
DM-DE coupling, and its discovery could mark a step forwards in the understanding of their 
nature. 
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Here we test this possibility by performing some numerical experiments. We assume 
DDE and CDE due to a SUGRA potential and use MCMC techniques to fit the following 
parameter set: ojo^ci ^o,b, t, ^, "ng, and (constant) w; 9 is the angular size of the sound 
horizon at recombination (see however below), Ug and ^4^ are spectral index and amplitude 
of scalar waves, no tensor mode is considered. 

The plan of the paper is as follows: In Section 2 we discuss how artificial data are 
built, outlining the models selected, the DDE potential used and the sensitivity assumed. 
In Section 3 we briefly debate the features of the MCMC algorithm used and illustrate a 
test on its efficiency, also outlining the physical reasons why some variables arc more or less 
efficiently recovered. In Section 4 we discuss the results of an analysis of DDE artificial data, 
against the w — const, assumption. In Section 5 we briefly summarize why and how CDE 
models are built and do the same of Sec. 4 for CDE models. This section yields the most 
signiflcant results of this work. In Section 6 we draw our conclusions. 



2. Building cirtificial data 

In order to produce artificial data we use CAMB, or a suitable extension of it (sec 
below), to derive the angular spectra C[^'^\ C\^^\ C\^^^ for various sets of parameter 
values. Artificial data are then worked out from spectra, according to Perotto et al. (2006) . 

The analysis we report is however mostly based on a single choice of parameters, 
WMAP5 inspired (see Komatsu et al., 2008, Spergel et al. 2007), in association with three 
different /5's. Two other parameter choices will also be used: (i) to test the efficiency of 
the MCMC algorithm; (ii) to add results for a still lower (3 value, so strengthening our 
conclusions. 

The parameter sets for the most general case and the (ii) case are shown in Table 1. In 
these cases DE is due to a scalar field self-interacting through a SUGRA potential 

V{<t^) = (A"+Vr ) exp(47r 0VmJ) , (1) 

[rup : Planck mass) which has been shown to fit CMB data at least as well as ACDM 
(Colombo & Gervasi 2007); as outlined in the Table, we input the value of the energy scale 
A; the corresponding a value is determined by the program itself; the choice A = 1 GeV is 
consistent with Colombo & Gervasi (2007) findings and is however scarcely constrained by 
data. The meaning of the coupling constant j3 is discussed at the beginning of Section 5 . 

Starting from the spectral components and assuming that cosmic fluctuations are dis- 
tributed according to a Gaussian process, we generate realizations of the coefficients of the 
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Model A 


Model B 




1 


1 


10^ iUo,b 


2.273 


2.4 


^o,c 


0.1099 


0.11 


100 /i 


71.9 


85 


^opt 


0.087 


0.15 


ln(10i° As) 


3.1634 


3.1355 


Us 


0.963 


1 


A/GeV 


1 


1 









/? 


0.05 


0.02 




0.10 





Table 1: Cosmological parameters for artificial CMB data. 



spherical harmonic expansions for the temperature and E-polarization fields, according to 
the expressions 

CTslli . (2) 
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where both (7^^, for any / and m, are casual variables, distributed in a Gaussian way with 
null averages and unit variance, so that the equality (glllgim) = is approached when the 
number of realizations increases. Together with eqs. ([2]) and ([3]) this grants that 
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if averages are taken over an "infinite" set of sky realizations. iFiom a single realization, we 
can define estimators of the power as 
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Taken independently from each other, cf'"'^^ and cj^^^^ at a given / are the sum of the squares 
of 2/ + 1 Gaussian random variables, so they are distributed according to a with 2/ + 1, 
which approaches a Gaussian distribution centered around the fiducial Ci^'^\ C\^^\ values. 
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as / increases. If we consider T and E simultaneously, the set of estimators cf"^\ Ci^^\ 
Cf^^ follow a Wishart distribution (see, e.g., Percival & Brown 2006). 

We consider here an idealized full-sky experiment characterized by: (i) finite resolution, 
that we shall set through a Full Width Half Maximum angle ■&fwhm, where the antenna 
sensitivity is reduced to 50 %, assuming a circularly symmetric Gaussian beam profile; (ii) 
background noise due to the apparatus, that we shall assume to be white, i.e. /-independent 
with an assigned variance a'j^. 

More in detail: for what concerns the Gaussian and circularly symmetric beam profile, 
its window function reads 

^ g-[KW)-.>V2 ^ with a = dFWHM/VS]^, (6) 

s being the spin of the signal (0 for anisotropy or i?-polarization spectra; its value however 
matters just for the lowest Z's). For what concerns noise, we shall consider the coefficients 
film such that 

with Nf^^"* = aj^6^^ , and R and S stand for cither T or E. As the sum of two indepen- 
dent Gaussian random variables is still Gaussian distributed, the statistics of the (beam- 
convolved) GMB -|-white noise field are given by Cl^^^ — C^^^^ \oBi\'^+Ni^^\ In an analogous 
manner from Cl^^^ we can define C\^^\ In the following we will consider both fiducial, i.e. 
C[^^\ or realized, C[^^\ model data sets. 

Under these simplified assumptions, the characteristics of an experiment are completely 
defined by the values of '&fwhm and ajf. Here we shall take '&fwhm — 7.0', while erf. ~ 
3 X 10~^(//K)^ and Up ~ 6 x 10~^(/iK)^. These can be considered conservative estimates of 
sensitivity in the forthcoming Planck experiment. 



3. The MCMC algorithm 

When the number J\f of the parameters to be determined from a given data set is large, 
the whole jV'-dimensional parameter space cannot be fully explored within a reasonable 
computing time. MonteGarlo Markov Ghain (MGMG) algorithms are then used, whose 
efficiency and reliability have been widely tested. 

In this work we used the MGMG engine and statistical tools provided by the GosmoMG 
package (Lewis & Bridle 2002; http:// cosmologist.info/cosmomc) . This tool set allows us to 
work out a suitable set of M Markov Ghains and analyze their statistical properties, both in 
the full jV'-dimensional parameter space, and in lower dimensional subspaces. Of particular 
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interest are the marginalized distributions of each parameter, obtained by integrating over 
the distribution of the other A/" — 1 parameters. 

The above algorithms implement the following steps, (i) Let Ci be a point-model in 
the parameter space. The spectra Ci of such model are computed and their corresponding 
likelihood Ci is evaluated according to 

where dCi stands for either fiducial or realized spectra. The above expression holds for a 
single TT or EE spectrum; for its generalization to 3 spectra and a discussion of the effects 
of anisotropic noise, sky cuts, etc. see, e.g. Percival & Brown (2006). (ii) The algorithm then 
randomly selects a different point (2, in parameter space, according to a suitable selection 
function. The probability of accepting (2 is given by min(£2/'Ci, !)• If C2 is accepted, it 
is added to the chain, otherwise the multiplicity A^^i of Ci increases by 1. (iii) The whole 
procedure is then iterated starting either from (2 or from (i again, until a stopping condition 
is met. The resulting chain is defined by [NiQ), where Q are the points explored and the 
multiplicity Ni is the number of times Q was kept. 

The whole cycle is repeated until the chains reach a satisfactory mixing and a good 
convergence. The first requisite essentially amounts to efficiently exploring the parameter 
space, by quickly moving through its whole volume. In particular, the fact that the algorithm 
sometimes accepts points with lower likelihood than the current point, avoids permanence in 
local minima, while a careful choice of parametrization and of the selection function minimize 
the time spent exploring degenerate directions. A good convergence instead guarantees 
that the statistical properties of the chains (or of a suitably defined subset of the chains), 
correspond to those of the underlying distribution. Here we implement the convergence 
criterion of Gelman & Rubin (1992), based on i?-ratio computation: convergence and mixing 
are reached when i? — 1 -C 0.1 . Fulfilling such criterion requires N ~ 10^ points, in the cases 
considered here. 

The algorithm was preliminarily tested by using ACDM as true cosmology and the set of 
parameter values shown in Figure [TJ Such Figure reports the likelihood distributions found 
in this case, both for input parameters (with asterisk) and for several derived parameters as 
well. Besides confirming that input values are suitably recovered, mostly with fairly small 
errors. Fig. [T] shows a broad non-Gaussian distribution for w values, with little skewness but 
significant kurtosis. This is a well-known consequence of the geometrical degeneracy present 
in CMB spectra, and clearly shows the difficulty of CMB data to yield results on DE nature. 
The non-Gaussian behavior is even more accentuated for some derived parameter, as Qo,c, 
Ho and the Universe age . We shall recover analogous features in the next plots; here we 
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Fig. 1. — Marginalized a posteriori likelihood distributions (solid lines), when artificial data 
derive from a ACDM cosmology. Vertical lines yield the input parameter values. Derived 
parameter panels are marked by an asterisk. Dotted lines show the average likelihood dis- 
tributions. 
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want to stress that they do not derive from fitting a partially unsuitable parameter set. 

Figure [T] also shows how errors increase when passing to some secondary parameter, 
from primary parameters specifically devised to break degeneracies. For instance, the errors 
on 6 and uJo,c are ~ 0.03% and ~ 1.3%, respectively. From these parameters. Ho and f2o,c 
are derived, whose errors are ~ 15 %. 

The point here is that the primary variable is 

e = rs{a*)/D{a*) , (9) 

a* setting the peak of the quasi-Gaussian last scattering band (LSB). It is then easy to see 
that 



rJa 



da Cs(a) Ho^o', 



1/2 



so that 9 is apparently independent from Ho and fio,m as, in the former equation 

HHa) /,^^WM\ (n) 



HI r^o/m ^ a J \ a 



while, in the latter one 



JP{a) 



and flo,m only sets the normalization of the latter term at the r.h.s., so that g{l) = 1 . As 
a matter of fact, these equations exhibit a mild dependence of 9 on Uo^m and Uofi, as ao/aeq = 
2.41xl0Vm(To/2.726K)-4[1.681/(l+0.227iV,)] and (cjc)-^ = 3[l+(3/4)KV^„, J(a/ae,)(l+ 
0.227 A^i,)] (here To is the present CMB temperature and N^, is number of neutrino families 
in the radiative background). Information of Qo,m and, thence, on Ho can be obtained only if 
the factor g{ao/a), yielding the evolution of the ratio between DE and matter at low redshift, 
is under control. It is certainly so if a ACDM model is assumed; but, if we keep w as free 
parameter, the uncertainty on it, ranging around 60 % , reflect the difficulty to obtain the 
secondary parameters flo,m and Ho- 

More in general, eq. (fT2l) shows that CMB data are sensitive to a change of g{ao/a) 
causing a variation of LSB distance D*. On the contrary, 9 does not discriminate between 
different w{a) yielding the same D*. 

Incidentally, all that confirms that it is unnecessary to fix a* with high precision, or 
to discuss whether it is a suitable indicator of the LSB depth. Knowing a* with ~ 1 % 
approximation is quite suitable to this analysis. 
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4. Dynamical DE vs. w = const. 

Let us now discuss wliat liappens if DE state equation, in the "real" cosmology, cannot 
be safely approximated by a constant w . To build data, we use then a SUGRA cosmology 
and our own extension of CAMB, directly dealing with a SUGRA potential both in the 
absence and in the presence of DM-DE coupling. Artificial data worked out for Model A 
with /3 = were then fit to the same parameters as in Figure [H The whole findings are 
described in Table [2] (second column). 

In Figure [2] we add further information, comparing the value distributions in fits as- 
suming either a cosmology with w = const, or a SUGRA cosmology, then including A, as 
parameter, instead of w. The two fits agree well within l-o", among themselves and with 
input values. In the Figure we show the fiducial case. 

Let us briefiy comment the results in the absence of coupling, when the fitting parameter 
set includes w instead of A: (i) uJo,b and Uo^c are recovered without any bias, (ii) The same 
holds for 6, if a numerically refined value of a* is used; the expression of a*, as above outlined, 
has a precision ~ 1 %. Errors on a*, however, are negligible in comparison with errors arising 
from the scarce knowledge of w{a). 

As outlined in the previous Section, the main problem to fix flo,m, ^o,a, and Ho, resides 
in the difficulty to determine the state equation of the expansion source since DE becomes 
substantially sub-dominant or dominant. This also makes clear why the formal error in the 
determination of the above three parameters is wider when the fit assumes a constant w : 
varying w yields an immediate and strong effect on DE contributions and state equation, 
extending up to now; on the contrary, only huge displacements of the energy scale A cause 
significant variations to DE contributions and state equation within the family of SUGRA 
cosmologies. In other terms, when A varies, the part of the functional space that w{a) covers 
is not so wide. 

Accordingly, the greater error bars on flo,m, ^o,a, Ho, found when we use w as a param- 
eter, do not arise because we fit data to a "wrong" cosmology, but because varying w leads 
to a more effective spanning of the functional space of w{a). 

As a matter of fact, the reliability of the fit must be measured by comparing the like- 
lihood of the best-fit parameter sets. Likelihood values, given in the next Section, do not 
discriminate between the two fits. Otherwise, a direct insight into such reliability is ob- 
tained by looking at the width of the marginalized a posteriori likelihood distributions on 
primary parameter values. Figure [2] then confirms that such distributions, although slightly 
displaced, have similar width for all primary parameters. 



- 11 - 




Fig. 2. — Marginalized likelihood distributions on parameter values. Data built with an 
uncoupled SUGRA model. The parameter space for the fit includes either log(A/GeV) (solid 
lines) or a constant w (dashed lines). Vertical dotted lines show the input parameter values. 
Derived parameters plots are with asterisk. 



A general conclusion this discussion allows to draw is to beware from ever assuming 
that error estimates on secondary parameters, as ^o,m, ^o,b, and Ho, are safe. 

5. Coupled DE vs. w = const. 

In cosmologies including cold DM and DE, the equation obeyed by the stress-energy 
tensors of the dark components 

y(c) ^ ^ Q ^ (13) 

yields 

-riPc + Pde) = -HPc+Pc + Pde +Pde)- , (14) 
CLT a 

if their pressure and energy densities are pc, Pde and pc, pde] t is the conformal time and dots 
indicate differentiation in respect to r. Eqs. f|T3l) and f|T^ state that no force, apart gravity, 
acts between standard model particles and the dark components. The eq. fll3p is fulfilled if 
DM and DE separately satisfy the eqs. 

y(de) ^ CT^-)ct>^, , T^^) ^^.^^ = -CT^^^cf)^, , (15) 

where T^'^^^'^^^ are traces of the stress-energy tensors. C can be an arbitrary constant or, 
e.g., a function of itself; when C = 0, the two dark components are uncoupled and so are 
the cases we considered up to here. 

It is however known (Wetterich 1995, Amendola 2000, Amendola & Quercellini 2003) 
that self consistent theories can be built with 

'^?> rrip ^ 

so that (3 is used to parametrize the strength of DM-DE coupling. 

By using MCMC techniques, we explored cosmologies with (3 = 0.05 and 0.1 (Model 
A); we also report some results for (3 = 0.02 (Model B). We aim to show that, when j3 is 
excluded from the parameter budget, the values MCMC provide for some parameters can be 
significantly biased, even though (3 is small. 

In Table[2] (3rd-4th columns) we report the results of fitting the parameter set including 
constant w, in place of A and f3. The Table allows an easy comparison with the uncoupled 
case, when all parameters unrelated to DE are nicely recovered. On the contrary, already 
for j3 = 0.05, the input uJo,c value is formally 4.5 a's away from best-fit, in the fiducial 
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Fig. 3. — Marginalized likelihood distributions on parameter values. Data built with coupled 
SUGRA models with (5 — 0.05; /3 = 0.1 (upper, lower panels). The parameter space for fits 
includes log(A/GeV) and (5 (solid lines) or just a constant w (dashed lines). Vertical dotted 
lines show the input parameter values. Pay attention to the different abscissa units in the 
two panel sets. 
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Fig. 4. — Anisotropy spectra comparison. The 3 values of (3 considered are clearly distin- 
guished thanks to their different input normalization and to the line type indicated inside the 
upper frame. The dotted lines (hardly visible at low I's in the upper panel) are the spectra of 
the corresponding best-fit model when assuming no coupling and w = const. . In the lower 
panel the relative differences between input and best-fit models are compared with WMAP5 
error size. 
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Fig. 5. — The comparison made in the upper panel of the previous Figure is extended to ET 
and E-polarization spectra. Model differences are hardly visible at very low I and on the I 
values where the ET spectrum changes sign. 
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Fig. 6. — Luminosity distance vs. redshift in models with {3 =0, 0.05, 0.10 (solid lines). Each 
model is compared with the corresponding best-fit constants model (dashed line) and with 
a SCDM model (dotted line) with equal Q;, . 
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case, and more than 5 cr's in some realization. The effect is stronger for /3 = 0.1, yielding 
a discrepancy ~ 8 cr's. That coupling affects the detection of uJo,c is not casual: in coupled 
models a continuous exchange of energy between DM and DE occurs. 

Several parameters exhibit a non-Gaussian distribution, however. This is visible in 
Figure [3l where we also provide a comparison between likelihood distributions when (i) the 
fit includes just w or, instead, (ii) the parameter space includes A and (3. 

The comparison allows to appreciate that, in the (i) case, Uo^c is apparently much better 
determined than in the (ii) case. Appreciable discrepancies concern also the 9 parameter 
and they reflect onto the derived parameters f2o,b, ^o,m and ifo, whose distribution however 
appears significantly non-Gaussian. 

These results deserve to be accompanied by a comparison between the likelihood values 
in the different cosmologies. The following table shows Xefr = — 21n(£) for the fiducial cases 
only: 

constant w A & /9 



/3 = 0.175 0.887 

(3 = 0.05 0.214 0.332 

p = 0.10 0.778 0.251 
The likelihood values for realizations are systematically smaller, as expected, but confirm the 
lack of significance shown here. Given that coupled models have 1 additional parameter over 
the w = const model, differences Axis <~ 1 indicate that both models are an equally good 
fit of data. The conclusion is that spectral differences, between input and best-fit models, 
lay systematically below the error size. As a matter of fact, we are apparently meeting a 
case of degeneracy. 

This statistical observation is corroborated by a direct insight into angular spectra, 
provided by Figures H] and [51 We show the behaviors of the Ci spectra for the Model A, 
with P = 0, P = 0.05, P = 0.1, compared with the Ci for the best-fitting {w = const.) 
model. Differences are so small to be hardly visible. Figure H] refers just to anisotropy. In its 
lower panel the ACi/Ci ratio is compared with WMAP5 l-o" error size (AC; is the difference 
between the spectra of the input and best-fit models). Notice that shifts by 1 or 2 units 
along the abscissa, at large / values, would still badly cut-off the apparent ratios. Clearly, 
degeneracy can be removed only if the error size is reduced by a factor ~ 10 or more, in the 
region between the first deep and the second acoustic peak. 

A further example we wish to add concerns a still smaller coupling intensity, P = 0.02 
(Model B, Table [3]). At this coupling level, the MCMC meet all input parameter values 
within a couple of a's. A closer inspection of plots similar to Figure [3l not reported here. 
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shows that the probabihty of reahzations yielding ujo^c formally more than 3-cr's away from 
the true value, is still > 8 % : the genie is still not completely back inside the lamp. 

6. Other observables 

The degeneracy observed in CMB spectra could be broken off through different observ- 
ables. We plan to deepen this aspect in a forthcoming work, by building detailed data sets 
accounting from the dependence on cosmology of the expansion rate and growth factor. 

A first insight into the actual situations can be however gained through an inspection 
of Hubble diagrams and transfer functions. 

In Figure [6] we show the redshift dependence of the luminosity distances for models with 
/5 = 0, 0.05, 0.1 and compare them with the corresponding best-fit models with constant 
w, as well as with SCDM models with the same value of These plots clearly indicate 
that a fit with SNIa data would hardly allow any discriminatory signal: discrepancies from 
constant-ty model increase with /5; but, even for the /3 = 0.10 case, they hardly exceed 
~ 10 % of the difference from SCDM. 

In Figure [7] we then exhibit the transfer functions T(A;) (multiplied by /c^'^ to improve 
the visibility of details) for models A. 

In the cases /3 = and /5 = 0.05, we notice slight displacements for the BAO system 
and the slope. The actual setting of BAO's is however subject to non-linear effects and 
residual theoretical uncertainties are wider than the shifts in the plots. The change of slope 
is also easily compensated by a shift of by ~ 0.01, widely within expected observational 
errors. Notice then that the relative position of the input and best-fit functions is opposite 
in the two cases. Accordingly, in the intermediate case /3 = 0.033, not shown in the plot, 
the overlap is almost exact and the observed scale dependence of the growth factor, at low 
redshift, does not break the degeneracy at all. 

The situation is different when /3 > 0.05 is considered. For /3 = 0.1, the BAO displace- 
ment is indeed relevant; to compensate for the change of slope, we would then require a shift 
of Ug greater than 0.1 . 

A numerical analysis can therefore determine at which value, probably intermediate 
between 0.05 and 0.1, the CMB spectra degeneracy is broken. 
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7. Conclusions 

The nature of the dark cosmic components is unclear. DE could be a self-interacting 
field, yielding a scale dependent state parameter w{a) . Which bias does then arise on cos- 
mo logical parameter estimates, if performed by assuming w = const. ? This question regards 
also parameters which do not describe DE; their estimate could be biased because the true 
model is not directly explored. A first conclusion of this work is that such bias exists but, 
in the cases we treated, yields acceptable displacements, within l-cr. 

Suppose however that future data allow to exclude Ug — 1, within S-cr's, when we 
assume w — const. . Setting < 1 discriminates among infiationary potentials. It would be 
however legitimate to assess that, at that stage, such conclusions would still be premature. 

Our analysis was then extended to the case of DM-DE coupling. The idea that DM 
and DE have related origins or arise from the same field has been widely pursued (see, e.g. 
Kamenshchik, Moschella & Pasquier 2001, Bento, Bertolami Sen 2002, 2004; Mainini & 
Bonometto 2004; for a review, see Copeland, Sami & Tsujikawa 2006). Couphng causes 
DM-DE energy exchanges and this option was first explored in the attempt to ease the 
coincidence problem. 

Unfortunately, when a DM-DE coupling, strong enough to this aim, is added to models, 
predictions disagree with data. This does not forbid, however, that the physics of the dark 
components includes a weaker coupling or that a stronger coupling is compensated by other 
features (see, e.g., La Vacca et al. 2008). 

In this work we pointed out that a significant degeneracy exists, so that we can find an 
excellent fit of CMB data for coupled DE cosmologies just by using constant-w uncoupled 
cosmologies. The fit is so good that even likelihood estimates do not allow to distinguish 

between the "true" cosmology and the best-fit constant-w model, at the present or foresee- 
able sensitivity levels. Unfortunately, however, the values obtained for several parameters 
are then widely different from input ones. 

In fact, if we ignore the coupling degree of freedom, when data are analyzed, we can 
find biased values for some primary parameter as Uo^c, for which input values lay ~ S-ci's 
away from what is "detected". Also ^lo,m and Ho, which are secondary parameters, are 
significantly biased. In particular, with a coupling as low a.s (3 = 0.05, we found model 
realisations yielding Ho estimate ~ 2-0" 's away from the input value. If we keep to the 
fiducial case, however, the probability to find Ho > its input value is 36.5 % ior P — 0.05 
and 26.5% for /3 = 0.1. This outlines the tendency to find smaller values than "true" ones. 

As a matter of fact, however, because of the uncertainty induced by our ignorance on 
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DE state equation, the width of errors on secondary parameter significantly exceeds the 
width for primary ones, and a 5-a discrepancy is partially hidden by such ignorance. This 
agrees with the fact that data analysis shows a significant increase of errors on secondary 
parameters when the set of models inspected passes from ACDM to generic-tt; cosmologies. 
For instance, in the very WMAP5 analysis, the error on Hq increases by a huge factor ~ 6, 
when one abandons the safety of ACDM to explore generic constant-w models. A general 
warning is then that errors obtained assuming ACDM are to be taken with some reserve. 

An important question is whether the bias persists when different observables, besides 
CMB, are used. A preliminary inspection shows that SNIa data would hardly provide any 
discrimination. On the contrary an analysis of fluctuation growth can be discriminatory if 
P >~ 0.07-0.08 . 

An even more discriminatory signal could be found in the z dependence of the growth 
factor is suitably analyzed. Observational projects aiming at performing a tomography of 
weak lensing, like DUNE-EUCLID (see, e.g., Refregier et al. 2006, 2008), therefore, can be 
expected to reduce this degeneracy case. 
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Table 2: Results of an MCMC analysis, seeking the parameters listed in the first column, on 
artificial CMB data built with the parameter values also listed in the first column, but using 
a SUGRA cosmology, whose A and {3 are shown in the header. For each parameter, the first 
line yields results for the fiducial case, the next 3 lines for model realizations. 
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Fig. 7. — Transfer functions for the models A. The dotted lines yield the transfer function 
of the models with w— const., yielding the best-fit for CMB spectra. 
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Input liiodel: coupled SUGRA (A = 



Parameter & Input value 



In 



no,m = 0.185 



IGeV. J = 0.02) 
b.f. value ±a 



0.254- 0.057 + 0.058 
0.294 - 0.062 + 0.062 
0.268 - 0.059 + 0.059 
0.219- 0.054 + 0.054 



102a;o,6 = 2.400 



uJo,c = 0.1100 



T = 0.1500 



w {ect z — 0) 



Us = 1.0000 



(10^%) = 3.136 



100 /i = 85 



2.400 + 0.016 
2.424 + 0.016 
2.423 + 0.017 
2.391 + 0.016 
0.1127 + 0.0014 
0.1120 + 0.0013 
0.1121 + 0.0013 
0.1130 + 0.0014 
1.11135 + 0.00032 
1.11174 + 0.00031 
1.11170 + 0.00032 
1.11108 + 0.00032 
0.1515 + 0.0059 
0.1413 + 0.0067 
0.1412 + 0.0061 
0.1510 + 0.0061 
-0.70-0.16 + 0.16 
-0.65-0.17 + 0.17 
-0.67-0.17 + 0.16 
-0.63-0.14 + 0.13 
0.9998 + 0.0042 
0.9994 + 0.0039 
0.9990 + 0.0041 
0.9976 + 0.0041 

3.138 + 0.011 
3.119 + 0.013 
3.119 + 0.012 

3.139 + 0.012 
79-9 +10 
77 - 10 + 10 
77 - 10 + 10 
75-8 +8 



Table 3: As previous Table, with different input values of cosmological parameters and (3 as 
small as 0.02 



